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1. INTRODUCTION 


The coordinates of several tracking stations tied to the NAD datum were 
computed at OSU through available observations to the GEOS-I satellite. 

Up to date results of the NA6 adjustment [Mueller et al. , 1969] and NA8 
adjustment [Mueller and Reilly, 1971] had been presented. The latter 
solution was performed using height constraints deduced from the SAO-69 
geoid [Gaposchkin and Lambeck, 1970]. 

0 

Recently a new detailed geoidal map with claimed accuracies of ± 2 m 
(on land), based on gravimetric and satellite data, was presented f Vincent, 
Strange and Marsh, 1971]. With the new geoid and the orthometric heights 
given in [NASA, 1971] more reliable height constraints were calculated 
and applied. 

The basic purpose of this experiment was to compute the new solution NA9 
by defining the origin of the system, from the point of view of error propaga- 
tion, in the most favorable position applying "inner constraints" TBlaha, 1971], 
and imposing new weighted height constraints to all of the stations. The 
major differences with respect to formerly published adjustments can be 
seen in Table 1. 
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2. THEORETICAL BACKGROUND 
2. 1 Normal Equations for Optical Observations 

The set of optical observations from MOTS and PC-1000 data had been 
previously screened [Mueller et al. , 1969] and a set of reduced normal 
equations of the form 

NX + U = 0 

obtained. The same symmetric coefficient matrix was used. This is 
composed of 3x3 blocks of the form [Mueller, 1968] 

N** = EM *] 1 - EM *] 1 (EMjj ) 1 M *] 1 + P* 

3x3 j i 1 

N*i = -EM*j (EMj ] 1 )' 1 M tJ 

3x3 1 1 

where 

= By Pi] 1 BiV 

10 0 
0 -cos 6 0 

0 0-1 

S = Solar motion matrix = R 2 (-x p )R l (-yp)R 3 (Sidereal Time) 

and 

k, i denotes particular ground stations 

s denotes particular simultaneous event 

l denotes any ground station participating in an event 

L is the summation over all ground stations involved in event j 
i 

E is the summation over all events observed by ground station k and/or i 
j 

P* weight matrix associated with any particular ground station 

3x3 

Pu weight matrix of any observed direction 

3x3 
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Table 1 

General Information on the Geometric Adjustment 
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Finally, the vector of constant terms is expressed by 
4 = -Mj 1 r Xk° - (LMjj 1 y 1 LM u x x°] 

3 11 

where as usual, the superscript (°) denotes initial approximate values. 
2.2 Normal Equations for Range Observations 

The general form of the normal equations similarly to the optical case 
can be formulated as: 

NX + U = 0 

where the 3x3 blocks in N are now computed by T Mueller, 1968] 

Nkk ~ Pkj ~ EaJj Picj rra T 13 Pu a ij 1 1 ^3 Pkj akj 

3x3 J i 1 

i . N k ! = p kj (Ea T t j p tJ a lj ) 1 a 13 p tj a^] 

3x3 j i 

and the vector of constant terms having the form 
. U* = -£a&j Pkj v kJ 
where 

v kJ = Residual of any observed range from a particular station 
(resulting from a preliminary least squares adjustment 
of any simultaneous event with the stations held fixed) 

P u = Weight of any observed range r u 



and 

X°, Y°, Z° = Appropriate Cartesian coordinates in the average 
terrestrial system. 

Other subscripts and symbols have the same meaning as in the case of 
optical observations. 
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2. 3 Constraint’s Contributions to the Normal Equations 

Two alternative definitions exist for the term "constraints". The 
absolute constraints represent certain conditions which have to be ful- 
filled exactly and with no uncertainties and the relative constraints (or 
weighted constraints) which have the same characteristics as the observations. 
In general the contribution of the functional constraint equations 

G(X, Lc) - 0 

to the normal equations can be found bordering the normal equation matrix 
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from where after elimination of Kc n it is easy to find 

X-i + C„ Pc n C n ]X n + U,.! + C T n P Cn W n = 0 
or 

IX- ! + N° n ]x n + U n _ i + lf n = o (1) 

where ISr n and L f n are the contributions to the coefficient matrix and 
constant vector of the normal equation due to the application of constraints. 

The coefficient n - 1 represents the normal equations of the previous set 
(without constraints). 

After the constraints are added the normal equations will take the 
usual form 

N n X n + U n = 0 

and we are in the position to obtain the contribution from a new set of constraints. 

Constraints can be applied between two stations k and 1 or to a single station. 
The contribution of these constraints to the matrix N (3x3 blocks) and 
U (3x1 blocks) can be schematically expressed in two different ways. 

Assuming that the matrix P c is always diagonal it is possible to express: 
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a) Contribution to the normals due to the constraint applied tp station k 
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l£ = Ck Pc W c 


b) Contribution to the normals due to the constraint between stations k and 1 




c[ 


Ck 


Vfr 

















































i 













I 

_ 









V 


- 

SI 




N kk 


Nki 










i 













II' 






__ 




_ 



— 


_ 

H 


N- 


lk 


Nn 


It 


Uk C k T P c C k ; i C k P c C v 

l£ = C k T P c W*; l f x = Cl P c W" 
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These blocks obtained as indicated above for the corresponding case 
will be the only ones computed and added to the original normal equations 
as expressed by formula (1). 

2.31 Relative Position Constraints 

Relative position constraints were used in order to combine the normal 
equations from the optical and range mode. The relative constraints used 
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and their weights are given in Table 1. The expression for the combi- 
nation of normals can be written as follows. 


rN + N RC ]X + U + U RC = 0 
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where and computed from (2) are the contribution to the original 
normal equations (NX + U - 0). 

In this case the functional constraint equations are 


Xk - X t = AX 

Y k - Y x = AY 

Z k ~ = Az 


Therefore: 



^,RC 

= I ; 

C* 


3 x 3 

3 x 3 

3 x 3 

and 

U* 

= 0 ; 

U* 


3 x 1 


3 x 1 

where 
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Thus, the diagonal elements of Prc are added to each element of the 
diagonal of the blocks kk and 11 of the matrix of the combined normals 
N, and subtracted from the diagonal elements of the blocks kl and Ik 
of N. There is not contribution to the vector U. 

2. 32 Length Constraints 

The range observations introduced scale into our adjustment. In 
addition to the SECOR scale, a distance constraint was imposed between 
stations 3861 and 7043, whose updated Cape Canaveral datum coor- 
dinates were derived from the high precision geodimeter traverse in 
the eastern United States. Due to a recent correction in their coordinates 
a difference of 3 m from the previously used value (see Table 1) was taken 
into consideration ("Meade, 1972]. 

The functional constraint equation is 

G Ch (X, Lch) ~ 0 or 

r<x k - X,) 2 + (Y k - Y f + <Z k - Zj) 2 ]^ = I*i 



and L _ ct 2 i difi = Variance of the chord 

* C H — ' 2 

(To (jo - A priori variance of unit weight 

Then the contribution to the normals are obtained applying (2) 

Nkk = (C C k H ) T P CH Ck CH 

3 x3 

= (C C ! H ) t P ch C C ! H 

3x3 

N;" = (C C k H ) T Pch Ci H 

3 x3 
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Figure 



uf = (c*y p CH \^ H 

3x1 

if" = (Ck H ) T P CH W CH 

3x1 

The first three expressions in the above are added respectively to the 
blocks Njtt, N u and of N; the last two expressions are added respectively 
to the constant vectors U k and 

2. 33 Height Constraints 

If the height of the station k is to be constraint, then 
N k H k = (C k H ) T P H C k H 

3 x3 

where 

C k = [cosc/^ cosA^, cos^ sinA^ , sinc^] 

3x1 

and 



where and A^ are the approximate geodetic coordinates and Oh k 
is the variance of the height for station k (see Table 4). 

The constant vector U H can be computed from 

U k = (C k ) T P H W H where 
W H = K 

At all stations a weighted height constraint was imposed in order to 
provide a valuable strengthening to the geometric network. 

After the completion of the previous adjustments (NA6, NA8), a 
new geoid (see Figure 1) became available [Vincent, Strange and Marsh, 
1971]. With the interpolated geoidal undulations for each station and their 
orthometric heights given in TNASA, 1971] reliable height constraints 
were computed by an iterative process self-explained in Figure 2 and 
summarized below. 
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SAO - NA8 = (SAO-NAD) + (NAD-NA8) 


NA8 ■-* SAO = NA8 + (SAO - NA8) 
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Figure 2 (continued) 


From the initial values of the shifts SAO-NA8 (computed using the 
published shifts SAO-NAD and NAD-NA8 in TMueller and Reilly, 1971]), 
the initial NA8 rectangular coordinates were shifted to the SAO-69 
origin and geodetic coordinates (<p 0 , Xo , ho) computed. 

Ellipsoidal heights h* then were calculated for each station using 
the undulations (N) from [Vincent, Strange and Marsh, 1971] (except 
for stations 1033 Fairbanks, Alaska and 3404 Swan Island whose undu- 
lations were interpolated from the SAO-69 geoid), and the sea level 
heights from [ NASA, 1971]. With the original and X„, and with this 
new height h* a new set of rectangular coordinates (with origin nearer to 
the geocenter) was obtained (X*, Y* Z*). From this set of coordinates 
and that of the NA8 new shift parameters (AX t , AY t , AZ t ) were computed 
and used to transform the NA8 set into NA8i, which is nearer to the geocenter. 
Following this procedure iteratively, several shifts of this kind to the "geo- 
center" were performed until the sum of the undulation differences with respect 
to N approached zero. Through this process a "best" set of geocentric 
coordinates (GC) was obtained. Using these and the NAD coordinates of 
the stations, the shift parameters NAD-GC were calculated which when 
added to the GC coordinates resulted in Cartesian coordinates in the NAD 
system. From these the heights h, to be constrained in the NA9 solution 
were computed (see step 18 in Figure 2 and Table 4). 
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2.34 Inner Constraints 


Unlike in the NA8 adjustment, "inner constraints" (free adjustment 
technique) were also imposed in the NA9 solution in order to define the 
origin of the system in its most favorable position from the error propagation 
point of view fBlaha, 1971]. 

The functional inner constraints equations can be written as 
C l X = 0 

where 



iH 

1 

0 

0 

1 

0 

0 ! 

1 1. 
1 

0 

o 

1 

o 

M 

II 

0 

1 

0 

0 

■ 

1 

0 

1 

• I 0 

1 

0 


0 

0 

1 

0 

0 

1 1 

1 

1 0 

0 

1 


and C has as many 3x3 unit blocks as unknown points. X is the 
set of corrections of the approximate coordinates of the unknown points. 

As it is known the "inner constraints" condition yields after the adjust- 
ment 

Trace (L x ) = minimum 

where £ x denotes the variance -covariance matrix of the coordinates of the 
ground stations following that the "inner adjustment" is the pseudo inverse 
solution by definition. 

In the most general application when the "best" origin, orientation and scale are 
sought the matrix C 1 has the form 
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where as usual (°) expresses the approximate coordinates of the points. 

If we represent the normal equations with the contribution of all the 
constraints (except inner constraints) by 

[N + N* + + N“]X + U + 0 + if" + U* =0 or 

NX + U = 0 


then the inner adjustment can be obtained by bordering the coefficient 
matrix N of the normal equations as 

(cyi 


N 


X 

“K i 


-U 

0 


Upon the addition of any kind of constraint to the normal equations, it 
becomes necessary to consider also its contribution to LV T PV. The degrees 
of freedom change as well. In order to compute the proper variance of 
unit weight the latter must be taken into consideration. 
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3. THE SOLUTION 


With the specific constraints mentioned above the NA9 solution was computed 
using the general OSUGOP program (from step 19 down in Figure 2) [Reilly 
etal. , 1972]. 

The coordinates of the NA9 solution are presented in Table 2 with their 
corresponding standard deviations. The coordinates transformed to the NAD 
are shown in Table 3. 

Table 4 presents the constrained heights (h) used at each station, the final 
undulations (Nqc) and their differences (£N), compared with those published in 
[Vincent, Strange and Marsh, 1971]. (For the procedure see steps 23 to 26 in 
Figure 2). In column AN and in parenthesis, the differences published previously 
in [Vincent, Strange and Marsh, 1971] are also shown. It can be seen that only 
two stations show substantial disagreements (3903 Herndon, Virginia and 3407 
Trinidad). It seems likely that the orthometric height as given in r NASA, 1971] 
for station 3903 has a gross error. The discrepancy with respect to station 3407 
may be due to the fact that it is situated in the Caribbean, where large geoidal 
gradients being present, the interpolation is uncertain. 


4. COMPARISON WITH OTHER SOLUTIONS 

Table 5 summarizes the transformation parameters (systematic differences) 
between the NA9, the "GC", the NAD and the SAO systems (the latter from 
[Gaposchkin and Lambeck, 1970]). Two sets of parameters are listed. The 
first was obtained through the assumption that only translation exists between the 
sets of coordinates. The second set presents the general seven parameter solution. 
In this latter transformation, first the rotations were computed from direction 
cosines independent of the origins and scales. The scale factors were calculated 
from the chord distances independent of the origins and orientations. These 
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parameters (rotations and scale) constrained with their variances were used 
in the solution for all seven parameters [Kumar, 1972]. Tables 6 to 9 show 
the constrained solutions with the resulting variance-covariance matrices 
and the correlation coefficient matrices. In the variance-covariance matrices 
the angular units are in radians . 

5. CONCLUSIONS 

The solution for the station coordinates (Tables 2 and 3) appear to be 
completely valid. The standard deviations of the coordinates are all acceptable 
being more realistic than those of previous adjustments especially in the geodetic 
heights. The obtained origin for the "GC" system seems to be very near the 
geocenter, despite its geometric deviation. 
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Table 2 

Coordinates of the North American GEOS-I Tracking Stations from the NA-9 Free Geometric Adjustment 
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All coordinates and standard deviations in meters. 





Table 2 (continued) 
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Table 3 

NAD Coordinates of the North American GEOS-I Tracking Stations 
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Note: The above coordinates were arrived at by applying the shifts AX =--6. 7m, AY - 0.2m and 
AZ = 0.4m to the NA-9 coordinates and then converting these values to ellipsoidal coordinates on 
the ellipsoid a = 6378206. 4 m, b = 6356583. 8 m. 





Table 3 (continued) 
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Table 4 


Height Constraints and Undulations 
(all units in meters) 


Number 

- 

Constraints 

N 

AN 

Station 

h 

CT h 



1021 

Blossom Pt. , Md. 

9 

3 

- 27 

-26 

- 1 (- 7) 

1022 

Fort Myers, Florida 

23 

3 

- 16 

-18 

2 < 1) 

1030 

Goldstone, Calif. . 

898 

3 

- 23 

-27 

4 < 8) 

1032 

St. John's, Nswf. 

102 

5 

12 

13 

- 1 

1033‘ r 

Fairbanks, Alaska 

145 


16 



1034 

E. Grand Forks, Minn. 

251 

3 

- 13 

-18 

5 (11) 

1042 

Rosman, N.C. 

916 

3 

- 23 

-22 

- 1 (- 3) 

3106 

Antigua, W. I. 

12 

3 

- 45 

-40 

- 5 (- 2) 

3334 

Stoneville, Mississippi 

45 

3 

- 20 

-19 

- 1 

3400 

Colorado Springs, Col. 

2184 

5 

- 4 

-10 

6 

3401 

Bedford, Mass. 

89 

5 

- 27 

-21 

- 6 

3402 

Semmes, Alabama 

84 

3 

- 21 

-18 

- 3 (-12) 

3404 + 

Swan Island 

79 

7 

- 32 



3405 

Grand Turk, B. I. 

0 

5 

- 51 

-47 

- 4 (-27) 

3406 

Curacao, N. Antilles 

44 

5 

- 30 

-26 

- 4 

3407 

Trinidad, Tobago 

285 

5 

- 50 

-34 

-16 

3648 

Hunter AFB, Georgia 

19 

3 

- 19 

-24 

5 (- 5) 

3657 

Aberdeen, Maryland 

7 

3 

- 27 

-26 

- 1 (- 4) 

3861 

Homestead, Florida 

16 

3 

- 22 

-22 

0 

3902 

Cheyenne, Wyoming 

1882 

5 

- 8 

-10 

2 

3903 

Herndon, Virginia 

132 

3 

-100 

-26 

-74 

5001 

Herndon, Virginia 

132 

3 

- 27 

-26 

- 1 

5333 

Stoneville, Mississippi 

45 

3 

- 17 

-19 

2 

5649 

Hunter AFB, Georgia 

23 

5 

- 23 

-23 

0 

5861 

Homestead, Florida 

22 

3 

- 27 

-22 

- 5 (- 5) 

7036 

Edinburg, Texas 

72 

3 

- 8 

-11 

3 

7037 

Columbia, Missouri 

270 

3 

- 17 

-24 

7 ( 10) 

7039 

Bermuda 

26 

3 

- 38 

-36 

- 2 (- 2) 

7040 

San Juan, P. R. 

57 

5 

- 40 

-41 

1( 5) 

7043 

Greenbelt, Maryland 

56 

3 

- 30 

-26 

- 4 (-15) 

7045 

Denver, Colorado 

1787 

3 


-10 

3 ( 16) 

7072 

Jupiter, Florida 

26 

3 


-24 

1(1) 

7075 

Sudburg, Canada 

276 

3 


-31 

3 ( 20) 

7076 

Kingston, Jamaica 

473 

3 


-23 

3 ( 20) 


*The "geocentric" coordinates were obtained from the NA-9 by adding the 
following shifts: AX = -51. 7 m, AY = 144. lm, AZ = 210. 5 m. 

+ Undulations from SA0 69 geoid. 
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Table 5 


Transformation Parameters 



NA9-NAD 

NA9-"GC" 

NA9-SAO 

"GC"-NAD 


28 

34 

11 

32 

£ •««* 

AX(m) 

6.7 ± 1.3 

51.7 ± 0.7 

35.3 ±2.6 

-44. 8 ± 1. 3 

CJ -01 






U C 
C3 ci 

AY(m) 

-0.2 ± 1.1 

-144.1 ±0.8 

-148.3 ±2.6 

144.8 ±1.1 

(X u 






CO 

AZ(m) 

-0.4 ±1.2 

-210.5 ±0.9 

-175. 0 ± 2.6 

209.8 ±1.2 

* 

AX(m) 

-3.5 ±2.4 

40. 0 ± 1. 6 

25. 9 ± 6. 6 

-38.0 ±2.2 

w 

£ 

AY(m) 

-18.4 ±2.5 

-142.9 ±1.1 

-200.7 ±7.9 

125.5 ±2.3 

+-> 

AZ(m) 

19.1 ±2.5 

-204.1 ±1.5 

-172.5 ±7.6 

226.7 ±2.3 

f-i 

(D 

<U 

0 Z (") 

-0.83 ± 0. 06 

-0.34 ± 0.04 

-0.62 ± 0.19 

-0.36 ± 0. 06’ 

s 

cd 

9y( // ) 

0.53 ± 0.05 

-0.20 ± 0.03 

0.14 ±0.16 

0.83 ± 0.04 

G- 

9x( // ) 

-0.14 ± 0.08 

-0.21 ± 0.05 

0.94 ± 0.24 

-0.10 ± 0.07 


f(xl(f ) 

-4.23 ± 0.32 

-0.15 ± 0.05 

-6.90 ±1.20 | 

-4.04 ± 0.29 


♦Rotation and scale factor parameters constrained. 
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Table 6 

TRANSFORMATION PARAMETERS: NA9-NAD 
Scale Factor and Rotation Parameters Constrained to Independently Determined Values 

AX 

METERS 
-3.51 


0.5730+01 

0.1230+01 

0.1920+01 

-0.3740-07 

C .5260-06 

0. 2720-06 

-0 

0.1230+01 

0 • t>060+0 1 

0.7240+00 

0 .5120—06 

0.2220-06 

0.7580-07 

-0 

0.192D+01 

0.7240+00 

0 .6060+01 

— G . 36 4 D— 06 

0.2900-06 

0.9240-07 

-0 

— 0.374D-0? 

0.5120-06 

—0 .3640—06 

0.1000-12 

-0.8390-15 

0.383D— 17 

0 

0.5260-06 

0.2220-06 

0.290D-G6 

-0.8390-15 

C .9220—13 

0.1600-13 

-0 

0.272D-06 

0.7580-07 

0.9240-07 

0. 3830— 17 

C. 1600-1 3 

0.5410—13 

-0 

-0.3650-06 

-0.5060-06 

— 0 .7350—06 

0. 1030-14 

— 0.573D— 13 

-0.20 70-13 

0 


COEFFICIENTS OF CORRELATION 


0.1000+01 

0. 2080+00 

0.3260+00 

-C. 4930-01 

(.7230+00 

0.4890+00 

-0 

0.2080+00 

0.1000+01 

0.1190+00 

G.O57D+00 

0.2970+00 

0.1320+00 

— 0 

0.3260+00 

0. 119D+00 

0.1000+01 

-0.4670+00 

0.3P8D+00 

0.161D+00 

-0 

— 0 .4930— Cl 

0.6570+00 

-0.4670+00 

0.1000 +01 

—0 .8720—02 

0. 5210-04 

0 

0.723D+C0 

0.2970+00 

0.3680+00 

-0.6720-02 

0.1000+01 

0.2270+00 

-0 

0. 4890+00 

0.1220+00 

C .1610+00 

0.5210-04 

0.2270+00 

0.1000+01 

-0 

— C.400 D +00 

—0 . 5400+0C 

—0 .7850+00 

0.851 0-02 

-0 .4950+00 

—0 .2340+00 

0 


AZ 

€ 

Bz 

e y 

e„ 

METERS 

( 10.D-F5 ) 

SECONDS 

SECONDS 

SECONDS 

19.09 

—4.23 

-0.83 

0.53 

1 

c 

• 

VARIANCE 

-COVARIANCE MAI 

R 1 X 



3650— Ct> 
506D-06 
735D-06 
1030-14 
573D— 13 
2070-13 
1450-12 


.4C0D+00 
.5400+00 
.7850+00 
.8510-02 
.4950+00 
.2340+00 
. aOCD+OI 
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Table 7 


TRANSFORMATION PARAMETERS: NA9-GC 
Scale Factor and Rotation Parameters Constrained to Independently Determined Values 


AX AY AZ € 0j 0y 6* 

METERS MITERS METERS (10.P+-5) SECONDS SECONDS SECONDS 

40.00 -142.92 -204.13 -0.15 -0.34 -0.20 -0.21 


VARIANCE - COVARIANCE MATRIX 


0. 268 0 +01 

0.6210+00 

0 .6860+00 

—0. 1720—08 

0. 2240-06 

0 . 1920—06 

-O. 1470-06 

G.621D+00 

0.1290+01 

0 .9360 + 00 

0.1 38D— 07 

0.9630-07 

0.3840-07 

-0.1920— C6 

0 .6660 +00 

0 . 9360+00 

0. 2120+01 

-0. 8570-08 

0.1080-06 

G. 3470-0/ 

-0.2670-06 

— 0 .1720—1)8 

0.1380-07 

-C- .8570-08 

0 • 25fcl>— 14 

-0.470D-17 

0.6590-17 

0.1290-17 

0.2240-06 

0.9630-07 

G . 108D— 06 

-0.4700-17 

0.3760-13 

0. 6300-14 

-0.2160-13 

0.142D— 06 

0.3840-07 

0.397D-C7 

0.659D— 17 

0 .8300—14 

0.2950-13 

—0 . 1Q2D-13 

— 0.147D-06 

-0. 1920-06 

-0.2870-06 

0.129D-17 

— 0. 2160—13 

-0. 1020-13 

0.560D— 13 


COEFFICIENTS OF CORRELATION 


O.IOOD+Ol 

0.3350+00 

0.2380+00 

—0 • 2030— C 1 

0.7060+00 

0.505D+00 

-0. 3790+00 

0.335D+0G 

0.1000+01 

G .5670+00 

0.239D+00 

0 .4380+00 

0. 197D+0G 

-0.715D+00 

0.268D+00 

0.5670, +0 0 

0. 1000+01 

—0 • 1 160+00 

0.3830+00 

0 . 139D+00 

— 0.834D+00 

-0.2080-01 

0.2390+00 

-0.1160+00 

0 . 10CD+01 

-6 .9 770-03 

0 .7560—03 

0.107D— 03 

0.7060+00 

0 .4380+00 

0. 3830+00 

—0 .4770—03 

0. 1000+01 

G. 2490+00 

-0.4710+CG 

0.5050+00 

0.1970+00 

0. 1390+00 

0.756D-03 

C.249D+0G 

0.1000+01 

-0.2510+00 

-0.3790+00 

— 0 .7150+00 

-0.8340+00 

G.1C70-03 

-0.4710 +00 

— 0 .25 10+00 

0 . iOOD+u 1 
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Table S 


TRANSFORMATION PARAMETERS: NA9-SAO 
Scale Factor and Rotation Parameters Constrained to Independently Determined Values 

&x ay az < e* e y e* 

METERS METERS METERS < 10.0+5) SECONOS SECONDS SECONDS 
25.-39 -200 .72 -112 . 54 -6.40 -n.fc2 0.14 0.94 


VARIANCE - COVARIANCE MATRIX 


0.4910+02 

0 .800U+0 1 

0.1290+02 

— E' . 64 4D— 06 

0.5050-05 

0.3180-05 

-0.242D-05 

C*. 8000+01 

0.622D+0? 

-0.2310+01 

0 .7300—05 

0. 18 ID— 05 

0.618D-06 

—0. 4980— 05 

0.1290+02 

-9.2310+01 

0.5820+02 

-6.5330-05 

0.1P0U-05 

0.4220-Ob 

-0.6450-05 

—0.6440-06 

0.730 0—0 5 

—0.5 33b— 05 

\ 

0. 145D-1 1 

0.220D— 14 

0 .1440—14 

—0 .4320—14 

0.505D-05 

U. 1810-05 

G.130D-05 

0 . 2200-14 

0.F760-12 

0. I860— 12 

-0.3800-12 

0.3180-05 

0.618 i " — 0 6 

0.4220-06 

0. 1440-14 

0.186D-1? 

0.610D— 12 

-0.1430-12 

-0.2420-05 

— 0 .4980—0 5 

-0 .6450-05 

-0.4320-14 

-0. 3800-12 

-0.143D-12 

0. 131D-11 


COEFFICIENTS OF CORRELATION 


0.1 000 +01 

C. 1530+00 

0.2550+00 

—0 . 804D— 0 1 

0.811 D+CQ 

0. 6120+00 

-0.3180+00 

0.1 530+00 

0.1000+01 

—0.3 8 3D— 01 

C . 768D+00 

0 .2460+0 0 

0.100D+00 

-0. 5520+00 

0 .2550+00 

—0.363D— 01 

0.100U+01 

—0. 575D+00 

0 .2 52D+00 

0.708D-vl 

— 0.739D+0C 

— 0 .8040-01 

0.768D+0G 

— C.579D+C0 

•0. 1000+01 

C). 1950-02 

0. 1530-02 

— 0.3 13D— 02 

0.811Q+00 

0.2460+00 

0. 2520+00 

0. 1950-02 

O.lOCP+Oi 

0.254D+00 

— 0.354D+00 

0.6120+00 

0.1000+00 

0.7080-01 

0.153D-C2 

6 .254D+09 

0 . 100D+01 

-0.160D+00 

— ('■ . 3 1 8 0 +0 0 

-0.5520 +00 

-G.739D+00 

-0.313C-02 

— G .3 54D+00 

-0. 1600+C0 

0.100D+01 
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Table 9 


TRANSFORMATION PARAMETERS: "GC-NAD 
Scale Factor and Rotation Parameters Constrained to Independently Determined Values 


AX AY. AZ C 0i 6 y 0, 

METERS METERS METERS (10.D+5) SECONDS SECONDS SECONDS 

-3B.C3 1 2 5 . 6 8 22c. 72 -A. 04 -Kt.3& 0.B3 -C.10 


VARIANCE - COVARIANCE MATRIX 


0.6839+01 

0.9720+CG 

C. 1 5 3 0+0 1 

—0.36 3D— 07 

0 .6300— n 6 

0.2300-06 

— 0.288D— 06 

0 .9720+00 

0.5100*01 

C .6660*00 

0.6-660—06 

0.1 790-06 

0.6160-07 

-0.6100-06 

0.1530+01 

0.6fit>0+00 

0.5070*01 

— 0 . 3120—06 

0 .2270—06 

0.7180-07 

—0.5930—06 

— 0 .3630—07 

0. 4463— 0 e 

-0.3120-06 

0.06 70-1 3 

-0.5785-15 

0.2210-16 

0.6810-15 

G.430D-C-6 

0 . 1790—06 

0.2270-06 

-0.578D-15 

0.7500-13 

U. 1310-13 

-0.6690-13 

C. 2 300 -0 6 

0 .6160—0 7 

0.71bD— 07 

0.221D-16 

0.1310-13 

0.666D-13 

-0.1660-13 

-0.2880-06 

-D.610D-06 

—0 .5930—06 

0.6b 1 0—1 5 

— ">.6690-13 

-0.1660-13 

0.1 170-12 


COEEF IC1ENTS OF CORRELATION 


0.1000*01 

0. 1960+00 

0.3090+00 

-0. 5290-01 

9.7160*00 

0.4850+00 -0.3860+00 

0.1960+00 

o.ieoa+o i 

0 .9550-01 

0.6710+00 

0.2 900+00 

0.127D+00 —0.5310+00 

0.3090+00 

0.9550—01 

0 . 1000*0 1 

— 0 .6 710*00 

0.3690+00 

0.1680+00 -0.7710+00 

— 0 .5290 -0 i 

0.6710*00 

-0.6 710+00 

0. 1000+01 

—0.71 70—02 

0 .3680—03 

0.6770-02 

0. 7160*00 

0 .2 900*0 0 

6 . 3690+00 

—0 . 7170—02 

0.1000+01 

0.2210+00 -0.6600+00 

0 .685D+00 

0. 1270*00 

0 • 1 680+00 

0. 3680—03 

0.2210+00 

0.100D+01 —0.2260+00 

—0 . 3 8 6 D *0 0 

—0 .531 D +0 0 

-0.7710*00 

0.6770-02 

—O. 6809+00 

— 0 .226D+00 

0.1 G0D+G1 
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